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SUMMARY 

Through the use of theoretical predictions of fluid properties, and exper- 
imental heat transfer and thrust- measurements , the zones of laminar, transi- 
tional, and turbulent boundary-layer flow were defined for the NASA Lewis 
1030:1 area ratio rocket nozzle. Tests were performed on the nozzle at cham- 
ber pressures from 350 to 1000. psia. For these conditions the throat diameter 
Reynolds numbers varied from 3x1 0^ to 10x10^. The propellants used were gase- 
ous hydrogen and gaseous oxygen. Thrust measurements and nozzle outer-wall 
temperature measurements were taken during the 3-sec test runs. 

Comparison of experimental heat transfer and thrust data with the corre- 
sponding predictions from the Two-Dimensional Kinetics (TDK) nozzle analysis 
program indicated laminar flow in the nozzle at a throat diameter Reynolds 
number of 3.2x10)5 or chamber pressure of 360 psia. Comparison of experimental 
and predicted heat transfer data indicated transitional flow up to and includ- 
ing a chamber pressure of 1000 psia. Predicted values of the axisymmetric 
acceleration parameter within the convergent and divergent nozzle were consist- 
ent with the above results. Based upon an extrapolation of the heat transfer 
data and predicted distributions of the axisymmetric acceleration parameter, 
transitional flow was predicted up to a throat diameter Reynolds number of 
22xl0 5 or 2600-psia chamber pressure. Above 2600-psia chamber pressure, fully 
developed turbulent flow was predicted. 


INTRODUCTION 

The objective of this study was to determine the rocket engine chamber 
pressure at which the boundary-layer flow in the NASA Lewis 1030:1 nozzle can 
be considered to be laminar, and at which pressure it will be turbulent. This 
information is needed in interpreting the measured thrust performance of the 
nozzle and to accurately predict its performance for specific test conditions. 
In 1986, the performance of a 1030:1 rocket nozzle at a chamber pressure of 
350 psia was measured at the NASA Lewis Research Center (ref. 1). A subsequent 
study determined that the flow within the nozzle was laminar (ref. 2). There 
is a need to determine the chamber pressure which will produce fully developed 
turbulent boundary-layer flow in the nozzle since turbulent flow produces more 
viscous drag than laminar flow and occurs in most high Reynolds number nozzle 
flows. Thus, the zones of laminar, transitional, and turbulent boundary-layer 
flow need to be defined for this nozzle as a function of Reynolds number or 
chamber pressure. 

Tests were performed at chamber pressures from 350 to 1000 psia. These 
conditions correspond to throat diameter Reynolds numbers from 3x10^ to lOxlO^. 


This paper is declared a work of the U.S. Government and Is 
not subject to copyright protection In the United States. 



The propellants used were gaseous oxygen and gaseous hydrogen. Thrust measure- 
ments and nozzle outer-wall temperature measurements were taken during the 
3-sec test runs. The hardware tested was an optimally-designed nozzle with a 
25° half angle of convergence, a l-1n. throat diameter, a contraction area 
ratio of 4.223:1 and an exit area ratio of 1030:1. The copper throat section 
was water-cooled and extended to an area ratio of 30:1. The nozzle skirt sec- 
tion was made of 0.25-in. carbon steel and was heat sink. 

In this report, the experimental heat transfer and thrust measurements are 
first compared to those predicted for laminar and turbulent flow to determine 
the state of the boundary-layer development. Then, predicted flow conditions 
will be used in conjunction with results from previous relaminarization studies 
to estimate the boundary-layer development. The final conclusions will be 
based upon a combination of the experimental and predicted results. 


SYMBOL LIST 


A area 

c* characteristic exhaust velocity 

Cf skin friction coefficient 

Cp specific heat at constant pressure 

F thrust 

H shape factor, displacement thickness/momentum thickness 

K acceleration parameter 

k conductivity 

P pressure 

q" wall heat flux 

R local radius 

r radius 

Re Reynolds number 

T temperature 

t time 

u velocity 

x axial distance from the throat 

a diffusivlty 
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c 


nozzle .expansion area ratio 
n c * characteristic exhaust velocity efficiency, percent 

0 momentum thickness 

p viscosity 

v kinematic viscosity 

p density 

a standard deviation 

Subscripts 

ax ax 1 symmetric 

c chamber 

capsule test capsule 
cr critical 

e boundary-layer edge conditions 

exit nozzle exit 

1 Inner 

o outer 

vac vacuum 

w wal 1 

9 momentum thickness 

® freestream 

BACKGROUND 

Figure 1 illustrates boundary-layer development as it is assumed to occur 
along the contour of a small rocket engine thrust chamber. The flow Is assumed 
to be stagnant at the interface of the injector and combustor. From there, a 
laminar boundary layer develops which transitions to turbulent. This transi- 
tion may or may not be completed depending upon the combustor length. If the 
transition is complete, the turbulent flow may be suppressed by the effects of 
acceleration (relaminarization) . If there is complete suppression of turbu- 
lence, the boundary layer becomes laminar. As the flow continues through the 
nozzle, If the fluid properties produce a Reynolds number high enough to cause 
transition, it transitions to turbulent. 
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Transition of a boundary layer from laminar to turbulent flow occurs as 
the Reynolds number Increases. Determining the point at which transition 
begins depends upon the flow conditions and, even in this day of high technol- 
ogy, Is still considered a "black art." The key to determining transition Is 
to determine the stability of the boundary layer and, therefore, the conditions 
which cause stability or Instability. Overall, accelerated flows (dp/dx < 0, 
favorable pressure gradient) like those which occur In a rocket nozzle, are 
considered more stable than decelerated flows (dp/dx > 0, adverse pressure gra- 
dient) (ref. 3). Therefore, accelerated flows would tend to transition at a 
higher Reynolds number. Also, the transfer of heat from the fluid (gas) to 
the wall causes the limit of stability to be larger (ref. 3.). Therefore, 
based upon heat transfer, transition In a rocket nozzle would occur at a higher 
Reynolds number. Conversely, transition occurs at a lower Reynolds number on 
a rough wall than on a smooth wall. In general, the flow conditions In a 
smooth-walled rocket nozzle promote stability In the boundary layer. 

Relamlnarlzatlon Is one of the names given to the phenomena whereby 
boundary-layer turbulence is suppressed by flow acceleration (favorable pres- 
sure gradient). To show the relationship between boundary-layer development 
and acceleration, let the acceleration parameter K be defined as: 



Then, the momentum Integral equation for turbulent plane flows with constant 
density and no suction or blowing can be written (ref. 4): 


dRe A C f 

u-ak * r - K<1 ♦ H)Re e 


If K Is positive, and If the term with K In it Is larger than the Cf 
term, the turbulent boundary layer approaches an equilibrium or stable condi- 
tion. The value of Ree at the point of equilibrium depends on K where 
large values of K lead to low values of Reg. Therefore, at some value of 
K, the equilibrium value of Re© will be below the critical Reynolds number 
for transition at which point turbulence production ceases and a laminar bound- 
ary layer emerges. This value of K Is called critical K or K cr . For 
values of K > K cr , the boundary layer is laminar. For values of K < K cr 
and above zero, there may be some suppression of turbulence which produces a 
transitional-type boundary-layer flow. In this report, the term "relaminar- 
ized flow" will refer to flow In which there Is complete suppression of turbu- 
lence (laminar flow) and partial suppression of turbulence (transitional flow). 

The effects of relaminarization In rocket nozzles were first noticed 
around the 1 950 ' s as a reduction in nozzle wall heat transfer. During this 
period, very large efforts in rocket engine research were beginning. Essen- 
tial to these efforts was the measurement of combustion chamber and nozzle 
heat transfer rates, especially around the throat section where heat transfer 
rates are the highest. At that time, turbulent boundary-layer flow, which pro- 
duces roughly twice as much heat transfer as laminar flow, was assumed. How- 
ever, the measured heat fluxes were significantly lower than predicted in the 
throat section using this assumption. It was speculated that the reduction in 
heat transfer was due to the onset of reverse transition, or relaminarization 
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of the boundary layer. At that point, no other experimental evidence was 
available to support this theory. Therefore, several studies were conducted . 
to explore this phenomena and Its- effects. 

Morettl and Kays (ref. 5) measured heat transfer of low velocity air with 
accelerating turbulent flow along a flat plate with varying wall temperature 
and varying freestream velocity. In this study, they used the acceleration 
parameter K as a measure of the acceleration effect. The acceleration param- 
eter Is essentially the ratio of a geometric acceleration parameter to the 
local Reynolds number based on chamber diameter, and Is usually defined as: 


Morettl and Kays found that at moderate values of K, there was a lower aver- 
age turbulence Intensity than found without flow acceleration. In fact, at 
values of K greater than 3.3xl0~ 6 , the heat transfer decreased sharply and 
approached the value for a purely laminar flow. At K > 3 . 5x1 0 -6 , there was 
no further reduction In heat transfer. It Is Important to note that these 
results were for a heated wall where heat transfer occurred from the wall to 
the fluid. 

Back., Cuffel and Massler (ref. 6) performed boundary-layer and heat trans- 
fer measurements for air flow through cooled, conical nozzles. Velocity and 
temperature profiles across the boundary layer were determined by measuring 
total pressure and total temperature. They found that at the higher chamber 
pressures, the heat transfer matched that of a turbulent flow. At lower pres- 
sures, they found a reduction In heat transfer along the convergent section. 
From their results, they estimated K cr as 2xl0 -6 to 3xlO -6 where they defined 
K as: 


2 

P e u e dx 

The boundary-layer measurements showed a change In the velocity and temperature 
profile near the wall at large values of K. In fact, the profiles became 
laminar-like. With regard to rocket nozzle geometry, they found that there 
was a greater reduction in heat transfer in the throat section of a nozzle 
with a larger half angle of convergence, but the same contraction area ratio. 
Therefore, their results suggested that an increase in the convergent half 
angle increases the effect of acceleration. 

They also noted that there was a lag in the thermal response to accelera- 
tion. Their conclusions were that acceleration affects the velocity distribu- 
tion In a flow undergoing relaminarization more than it affects the temperature 
distribution. The implication, then, is that relaminarization effects occur at 
larger Reynolds numbers than would be Indicated by the measured heat transfer. 
Thus, although the measured heat transfer may be characteristic of a turbulent 
boundary layer, acceleration may be affecting the friction coefficient and 
hence the drag may be less than for turbulent flow. 

Schoenman (ref. 7) presented thermal data from four different thrust 
chambers with different propellants. The data from his tests and other small 
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engine tests suggested that below a throat diameter Reynolds number of 2x10^, 
all nozzles will have laminar boundary-layer flow. In the Reynolds number 
range between 2xl0 5 and 4xl0 5 , the boundary-layer characteristics may be lami- 
nar, transitional, or turbulent depending upon the geometry, surface roughness, 
and the combustion effects. 


One of the most significant studies was done by Nash-Webber (ref. 8). In 
this study, he flowed air over a flat plate and varied the imposed pressure 
gradient. In effect, he produced relaminarization in compressible, turbulent 
flow over an adiabatic wall in the ranges of Mach and Reynolds numbers typical 
of rocket nozzles. During these tests, wall shear stress was measured. From 
his results, he determined the coupled effects of acceleration and characteris- 
tic Reynolds number. It should be noted that Nash-Webber varied from previous 
studies and used adiabatic wall conditions to define his acceleration parameter 
and momentum thickness Reynolds number as follows: 


K = 


p w 

p w u e 


dx 


and 



p Q U 0 

r e e 
p w 


He presented his results in the form of a "laminarization map" which defined 
the regions of complete suppression of turbulence (laminar flow), partial sup- 
pression (transitional flow), and no suppression (turbulent flow) for boundary- 
layer flow over an adiabatic wall. The purpose of the map was to plot the 
predicted boundary-layer conditions along a given nozzle wall as a trajectory 
on this diagram. Entry of the trajectory into the region defined for relami- 
narization indicates the onset of laminarization for that device. 


Another significant study (ref. 9) was done in which heat transfer and 
boundary-layer measurements were obtained for hot air flow through the conver- 
gent section of a cooled nozzle with a 30° half angle of convergence. The 
study was performed over a range of chamber pressures corresponding to throat 
diameter Reynolds numbers of 6x1 0 5 to 50xl0 5 . Boldman found that the throat 
heat transfer at a throat diameter Reynolds numbers below about 10 6 was equal 
to that predicted for a laminar boundary layer. He also found that at the 
lower chamber pressures, although the boundary-layer profile upstream of the 
convergent section indicated turbulent flow, a profile taken in the convergent 
section was altered suggesting transition towards a laminar-type boundary 
layer. Within this study, an axisymmetric acceleration parameter was derived 
from the integral momentum equation for axisymmetric flow. This parameter 
K ax was defined as 


v du 

f 00 oo 

' ax = 7dT 
00 


v dr 

+ 0.352 — 

u r dx 

co 


A critical value of K ax , Kax cr » was estimated using the criteria for flat 
plate boundary-layer transition and was found to be 2.88xlO -6 . Boldman found 
that calculations of K ax could be used to predict the depressions in heat 
transfer due to flow acceleration and relaminarization. However, it was noted 
that although K ax was less than K aXcr at the throat for the lower chamber 

pressures, laminar heat flux was measured. Therefore, it appears that high 
values of K ax upstream affected the boundary layer downstream to the point 
of lowering K ax . 
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Following the above mentioned studies and several others, it was deter- 
mined that for rocket nozzles relaminarization is controlled by the chamber 
contour and by the product of chamber pressure and thrust at the throat plane, 
since this product determines the throat diameter Reynolds number (ref. 10). 
Large convergence angles, large contraction ratios, and small throat radii of 
curvature promote relaminarization. It was also found that high values of K 
upstream of the throat could lower K cr near the throat. 

A summary of results from several studies appears in reference 10. A sug- 
gested upper limit for relaminarization Is set at a value of chamber pressure 
times thrust less than 10^ lbf 2 /1n. 2 . Further, the transition zone, or throat 
diameter Reynolds number range where the relaminarization process is started 
but not completed, is estimated to extend from 2x1 0 5 to 1 7x1 0 5 , depending on 
the convergence angle. 

As mentioned previously, the significance of studies on boundary-layer 
development and relaminarization is to enable rocket engine designers to pre- 
dict heat transfer in the combustion chamber and throat section. This informa- 
tion is necessary to design engine cooling systems. Currently, the concept of 
a high performance, compact (high chamber pressure) space engine is being eval- 
uated. This concept Includes the use of large area ratio nozzles. Because of 
the large area ratio, boundary-layer drag is a significant factor in evaluat- 
ing the design for optimal performance. According to prediction, turbulent 
boundary-layer drag reduces the thrust performance for the NASA Lewis 1030:1 
nozzle about 3.5 percent below the performance predicted for a completely lami- 
nar boundary layer (ref. 11). Therefore, it is important to understand the 
boundary-layer development in large nozzles such that the correct drag and, 
thus, correct performance can be predicted for specific run conditions. 

Tests have been performed on a 1030:1 nozzle contour at a low chamber 
pressure which resulted in laminar flow as evidenced by the experimental nozzle 
wall heat transfer. Tests need to be done on the same nozzle at a pressure 
which produces fully developed turbulent flow. To do these tests, it is neces- 
sary to define the chamber pressure at which acceleration has no effect on the 
boundary layer in the nozzle. Not only will this help to define an experimen- 
tal program, but it will also provide insight into the design of similar space 
engines. 


APPARATUS AND INSTRUMENTATION 
Test Facility 

Testing was done in the altitude test capsule at the NASA Lewis Rocket 
Engine Test Facility (RETF). Figure 2 is an illustration of RETF with cutaway 
views of the test capsule and spray cooler. The operation of the facility was 
as follows: When the engine was fired, the exhaust gases flowed into the dif- 

fuser where the kinetic energy of the exhaust was used to accomplish some of 
the altitude pumping. From the diffuser, the exhaust gases flowed into the 
spray cooler, where approximately one-half of the exhaust was condensed to 
water and flowed out the vertical drain line. The other half was pumped by 
the gaseous nitrogen ejectors shown mounted on top of the spray cooler. The 
exhaust from these ejectors was directed up through two short stacks and vented 
to the atmosphere. 
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The thrust stand used in this facility was capable of measuring 3000 lb. 
full scale and was attached to a foundation that was separate from the test 
capsule bulkhead. The thrust stand was designed to have a 2 a (standard devia- 
tion) variation of less than +0.1 percent of full scale, and was calibrated 
against a load cell that had a 2 a variation of less than +0.05 percent of 
full scale. 


Test Hardware 

The NASA Lewis 1030:1 thrust chamber hardware consists of three compo- 
nents: the combustion chamber, throat section, and nozzle skirt. Figure 3 is 

a sketch of this hardware. As shown, the combustion chamber was 6.0 in. long 
and had an internal diameter of 2.055 in. This combustor was made of copper 
with 1 .0-in. -thick walls and was uncooled or heat-sink hardware. 

The throat section had a half angle of convergence of 25° and a contrac- 
tion ratio of 4.223. The radius of curvature upstream was 1.0 in. and down- 
stream was 0.2 in. The throat section converged to a 1.00 in. throat and then 
diverged to an area ratio of 29.9. For the 350-psia chamber pressure tests, 
the throat section was made of nickel with a ceramic coating on the inner wall 
and a water jacket around the throat that extended to an area ratio of approxi- 
mately five. For all other tests, the throat section was made of copper with 
water coolant channels along the entire length of the section. This difference 
in hardware between the test cases should not affect the results. The 350-psia 
chamber pressure cases were found to have laminar flow in the nozzle and, 
therefore, the use of the water-cooled copper throat during these tests would 
not have altered the results. 

The nozzle skirt section extended from an area ratio of 29.9 to 1030. The 
axial length was approximately 4 ft and the exit plane diameter was 32 in. 

This section was made of 0.25-in. thick carbon steel and was uncooled or heat- 
sink hardware. The roughness measurement for the skirt was found to be 32 rms, 
which is in the regime of smooth surfaces. 

The nozzle contour past the throat was designed using the Rao Nozzle Con- 
tour Program (ref. 12) and the Boundary Layer Integral Matrix Procedure 
(BLIMP-J) Program (ref. 13). The Rao program calculated the inviscid optimal 
nozzle contour for a GH2/G02 1000:1 nozzle with a l-1n. throat, a chamber pres- 
sure of 1000 psia, and a nozzle length equivalent to 100 percent of a 15° cone. 
Using the Rao contour, the BLIMP-J program was run using GH2/G02 propellants, 
a chamber pressure of 1000 psia, and a mixture ratio of six, to calculate the 
boundary-layer displacement thickness along the nozzle. The BLIMP-J program 
then added this thickness to the original Rao contour to produce the nozzle 
coordinates that were used to make the hardware. 

Figure 4 is a photograph of the inside of the test capsule with the 1030:1 
nozzle in the process of being installed. 


Instrumentation 

Nozzle outer-wall temperatures were measured by chromel-constantan and 
chroma! -alumel thermocouples spot-welded to the nozzle surface. These thermo- 
couples were referenced to a 150 °F oven. A calibration point was taken on 
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the thermocouples prior to testing. The results Indicated all the thermocou- 
ples were within +3 °F accuracy. Four thermocouples were equally spaced cir- 
cumferentially at axial locations corresponding to area ratios of 50, 100, 200, 
300, 388, 500, 635, 800, and 975. Thermocouple outer wall measurements were 
not taken on the throat section since It was cooled hardware. 

Pressure taps were placed In the nozzle wall at the same axial positions 
as the thermocouples. These static pressure measurements were used to affirm 
that the flow In the nozzle was fully expanded. 


ANALYSIS 

Fluid properties, nozzle-wall heat transfer, and thrust performance were 
predicted using version 2.4 (December 1984) of the Two-Dimensional Kinetics 
(TDK) nozzle analysis computer code (ref. 14). TDK calculates the two- 
dimensional and kinetic effects on the performance of liquid propellant exhaust 
nozzles. The basic method of analysis used is the method of characteristics 
for kinetic expansions. A Boundary-Layer Module (BLM) has also been included 
In this program which takes the inviscid fluid properties predicted with the 
method of characteristics and calculates wall-heat flux and boundary-layer 
drag. Within BLM, the continuity, momentum, and energy equations for compres- 
sible boundary-layer flow in a symmetric nozzle are solved using a two-point 
finite difference method developed by Keller and Cebecl (ref. 15). The turbu- 
lence model used is the Cebeci-Smith eddy-viscosity formulation (ref. 16) which 
has been tested for many types of flows with various boundary conditions. 

Hardware specifications and test conditions were used to write the input 
files to TDK so that this program would accurately model the nozzle perform- 
ance. The experimental results used in the input files were: chamber pres- 

sure, mixture ratio, fuel injection temperature, oxidizer injection temperature 
and experimentally determined nozzle inner wall temperatures. 

Boundary-layer edge conditions in the combustor were input to the program 
based upon estimated values since no data was taken in the combustor. Wall 
temperatures in the throat section were also input based upon estimated values. 
To determine the sensitivity of the estimated wall temperatures in the throat 
section on the predicted results, the TDK program was first run with the esti- 
mated wall temperatures. Then, it was run again with a 25 percent increase in 
the wall temperatures. It was run a third time with a 25 percent decrease in 
wall temperature under the original estimate. The results showed that the pre- 
dicted acceleration parameter was not affected by the wall temperature, since 
it is calculated using only freestream properties. The predicted momentum 
thickness Reynolds number was affected by wall temperature to the extent that a 
25 percent change in wall temperature produced a maximum change of 2.7 percent 
in the momentum thickness Reynolds number along the throat section. A 25 per- 
cent change in throat-wall temperature produced a 0.015 percent change in pre- 
dicted thrust. Neither of these effects are significant for this study. 


DATA REDUCTION 

During each 3-sec test run, thrust measurements, nozzle wall static pres- 
sure, and nozzle outer-wall temperatures were taken. These measurements were 
taken at a rate of 50/sec, averaged in groups of five, and displayed at 0.1 sec 
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Intervals. The nozzle wall static pressure measurements Indicated that the 
nozzle flow was fully expanded In about 2.5 sec. All the data presented In 
this report was collected Immediately before the engine was shut down. All the 
data for the 350-psla chamber pressure tests are from references 1 and 2. 

Thrust was measured using a specially designed thrust stand which could 
measure 3000 lb. full scale and had a 2 a variation of less than +0.1 percent 
of full scale. To account for a zero shift due to variations In ambient (cap- 
sule) pressure, a correction was added to the measured thrust as a function of 
the capsule pressure. Then, vacuum thrust was calculated as: 

Fvac = F + (Pcapsule * A exit^ 

As mentioned previously, nozzle outer-wall temperatures were measured at 
several axial locations. Heat flux and Inner-wall temperatures at each thermo- 
couple location were determined by solving the energy equation: 

"S (&) ■ 7 aF [ <kr> If 

Assuming constant properties, no axial or circumferential conduction, adiabatic 
outer-wall and a linear rate of temperature rise, the following equations were 
obtained: 



Thus, Inner-wall temperatures and heat flux to the wall were calculated using 
the measured outer-wall temperature T 0 , the measured temperature rise rate 
(9T/8t), and properties of the carbon steel wall. 

As mentioned previously, four temperature measurements were taken at each 
axial location. For each of these measurements, a heat flux was calculated. 
These four values were then averaged to determine the heat flux at that axial 
location. Also, in order to account for the effect of combustion efficiency, 
the fluxes were adjusted in the following manner: 



The term (n c *)2 represents the adjustment for combustion efficiency based on 

the fact that c* is proportional to the square root of the total chamber 
temperature. 
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RESULTS 

The purpose of this study is to define the zones of laminar, transi- 
tional, and turbulent boundary-layer flow for the NASA Lewis 1030:1 rocket noz- 
zle. These zones will be established in terms of the throat diameter Reynolds 
number, which is determined by the chamber pressure and throat size. The pro- 
cedure used in this paper to determine the flow zones is as follows: First, 

the rocket nozzle was tested at various chamber pressures whi le measuring 
thrust and nozzle outer wall temperature. Static pressure measurements were 
also taken and indicated ful ly expanded flow for the test runs discussed in 
this report. Then, the TDK program was run to match each test case assuming 
(1) laminar flow which transitions to turbulent in the combustor, and (2) com- 
pletely laminar flow throughout. Next, the experimental wall heat flux was 
compared to the predicted fluxes for the two theoretical cases to infer the 
state of the boundary- layer development. Also, the thrust was compared to pre- 
dictions to determine if the drag produced compared with that of a turbulent 
or laminar boundary layer. Lastly, predicted fluid properties were used in 
conjunction with results from previous relaminarization studies to estimate 
the state of the boundary layer. 


Heat Transfer 

Figures 5 to 9 are plots of the experimental and predicted wall heat 
fluxes along the nozzle skirt for similar mixture ratios and chamber pressures 
of 360, 656, 747, 984, and 1004 psia, respectively. As shown, at 360-psia 
chamber pressure, the measured fluxes fall on the plot of the TDK prediction 
for laminar flow. As the chamber pressure increases, the experimental fluxes 
approach the turbulent plot. These results infer that at a chamber pressure 
of 360 psia, the flow was laminar. As the chamber pressure increased, the flow 
became transitional. In figure 9, at a chamber pressure of 1004 psia, the 
experimental fluxes fall about one-half way between the turbulent and laminar 
plots until near the end of the nozzle where the fluxes approach the turbulent 
plot. This Infers that the flow was still transitional, but was almost fully 
developed turbulent flow when it reached the exit plane. 

As a cross-check, figures 10 to 13 were made. The symbols represent the 
experimental heat fluxes for similar mixture ratio tests and chamber pressures 
of 354, 360,- 656, 747, 984, and 1004 psia. Each figure was plotted for one 
area ratio or axial location along the divergent nozzle. The dashed line 
represents the best fit to the data. 

In figures 10 to 13, a power curve fit was applied to the data over the 
range of chamber pressures from the 350 to 1004 psia and the following results 
were obtained: 


n 

Power curve equation: q a (P c > 

Area ratio n 

100 0.96 

300 .97 

•' 388 .91 

800 .83* 

*From a chamber pressure of 
656 to 1004 psia. 
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For external flows over a flat plate, the heat flux from a laminar boundary 
layer is proportional to the Reynolds number to the 0.5 power. This means that 
laminar boundary- layer heat flux should vary as a function of the square root 
of the chamber pressure. Similarly, for turbulent flow over a flat plate, the 
heat flux is proportional to the Reynolds number, or chamber pressure, to the 
0.8 power. As shown, at an area ratio of 100 and 300, the heat flux varies 
almost linearly (exponent is approximately 1) with the chamber pressure. This 
indicates that over this chamber pressure range, the heat flux Increased at a 
greater rate than expected with laminar or turbulent flow. This is the trend 
or slope one would expect to see with flow that is transitioning from laminar 
to turbulent. As the area ratio increases to 388, the exponent of the chamber 
pressure decreases indicating a decrease in slope of the data. However, the 
slope is still greater than that predicted for turbulent flow. At an area 
ratio of 800, the slope of the data between a chamber pressure of 656 and 
1004 psia is closed to that predicted for fully developed turbulent flow. This 
indicates that at the higher chamber pressures, as the flow approached the exit 
plane. It was almost fully developed turbulent flow. These results are con- 
sistent to those presented in figures 5 to 9. 


Thrust 

Figure 14 is a plot of the measured and predicted vacuum thrust for simi- 
lar mixture ratios. The upper line represents the predicted vacuum thrust for 
completely laminar boundary-layer flow. The lower line represents the pre- 
dicted vacuum thrust for boundary-layer flow that transitions to turbulent 
within the combustor. The symbols are thrust measurements taken at various 
chamber pressures. The uncertainty of the thrust measurements was calculated 
to be 1.12 percent (ref. 17). The predicted difference in vacuum thrust 
between laminar and turbulent flow is 2.9 percent. 

As shown, all the experimental data fall on or near the laminar predic- 
tion, indicating laminar flow. Based upon a one-to-one comparison of the data 
to the predictions, the following Information is obtained. At a chamber pres- 
sure of 360 psia, the thrust measurements are 2.2 percent above the laminar 
prediction and 5.1 percent above the turbulent prediction. At a chamber pres- 
sure of 656 psia, the thrust measurement is 0.3 percent above the laminar pre- 
diction and 3 percent above the turbulent prediction. At a chamber pressure of 
747 psia, the thrust measurement is 0.7 percent above the laminar prediction 
and 3.6 percent above the turbulent prediction. Thus, as the chamber pressure 
increases from 360 to 656 and 747 psia, the thrust measurements are lower rela- 
tive to the laminar and turbulent predictions. These results indicate an 
increase in drag relative to the assumed laminar flow found at 360 psia. This 
suggests the onset of transitional flow above 360-psia chamber pressure. 


Extrapolation of Experimental and Predicted Results 

Based upon the heat flux and thrust results, fully developed turbulent 
flow was not reached during the experiment. Therefore, the end of the transi- 
tion zone must be predicted based upon an extrapolation of the experimental 
data and predictions of flow conditions at higher chamber pressures. Figure 15 
is a plot of the heat flux obtained at an area ratio of 100 for similar mixture 
ratios. This point on the nozzle is about 7 in. downstream of the throat. The 
flux at an area ratio of 100 was chosen because it is the closest measurement 
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near the throat which is not affected by the water-cooled throat hardware. The 
heat flux is presented as a function of throat diameter Reynolds number. The 
upper curve represents the predicted flux for a turbulent boundary layer. The 
zone around the turbulent curve represents the approximate accuracy of the pre- 
diction (+20 percent). The lower curve represents the predicted heat flux for 
a laminar boundary layer. The experimental fluxes are represented by symbols 
with a least squares curve fit applied to the data. As shown, an extrapolation 
of the curve fit intersects the turbulent plot at a throat diameter Reynolds 
number of about 22x1 0 5 which corresponds to a chamber pressure of 2600 psia. 
Based upon the accuracy of the turbulent flux prediction, this' value is proba- 
bly accurate to within 500 or 600 psia. 

In order to predict the upper limit of the transition zone, the flow con- 
ditions obtained from the TDK program were used along with the calculated 
inner-wall temperatures to predict the boundary- layer development based upon 
the Nash-Webber laminarization map. The Nash-Webber map was chosen for. compar- 
ison because it defines the boundaries for laminar, transitional, and turbulent 
flow based upon the effect of acceleration. The values of K and Re© were 
calculated as defined by Nash-Webber for the 747-psia chamber pressure case. 

The distribution of K and Re© on the Nash-Webber map indicated that the 
flow was expected to be entirely turbulent. However, this result is not sup- 
ported by the heat flux and thrust data. In a previous study, Praharaj 
(ref. 18) also found that the Nash-Webber criteria predicted turbulence when 
the measured data indicated relaminarization. The problem appears to be that 
Nash-Webber' s study was based upon an adiabatic wall with his K parameter 
incorporating conditions at the wall. As mentioned previously, heat flux from 
the fluid to the wall causes the boundary layer to be more stable than heat 
flux from the wall to the fluid. Therefore, for a cooled wall or heat-sink 
hardware, the Nash-Webber laminarization limits would have to be lowered. 
Because the 1030:1 nozzle is heat-sink hardware, Nash-Webber' s criteria cannot 
be applied. 

Next, the predicted flow conditions were used to estimate the boundary- 
layer development based upon Boldman's results. Boldman's study was chosen 
because it was for axi symmetric flow through a cooled nozzle with a conver- 
gence half angle of 30°. These conditions are similar to those for the NASA 
Lewis 1030:1 nozzle. The acceleration parameter was calculated as specified 
in Boldman's paper (ref. 9). Local values of K ax for the nozzle are plotted 
in figures 16 to 20 as a function of the axial distance from the throat. 

The critical value of K ax used by Boldman was 2.88x10"^ and was based 
upon the assumption that laminarization occurs at a critical momentum thick- 
ness Reynolds number of 360, which is the classical value for transition on a 
flat plate. Later, Cebeci and Smith suggested that for compressible flow with 
a pressure gradient (as found with nozzle flow), the criteria for transition 
should be set at a Re© of 400 (ref. 16). Therefore, the critical value of 
K ax was recalculated using Boldman's derivation and a critical momentum thick- 
ness Reynolds number of 400. The critical K ax then becomes 2.5xl0 -6 . This 
is the value which is Indicated on figures 16 to 20. 

As shown, the value of K ax is quite high in the convergent part of the 
nozzle. However, as you cross the throat plane, the value drops. For each 
case, from a chamber pressure of 350 psia up to 1000 psia, one would expect 
laminar flow in the convergent section based upon K ax > K aX(;r , and laminar 

flow in the throat section based upon the upstream influence of K ax . This 
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matches Boldman's prediction of laminar flow In the convergent section below a 
throat diameter Reynolds number of 10 6 . The affect of acceleration when 
K ax < Kax cr Is predicted to decrease monotonical ly with the value of the 

acceleration parameter (ref. 8). Thus, in the divergent nozzle, one would 
expect less and less of an acceleration effect on the boundary-layer develop- 
ment as the flow expands to the exit. Therefore, one would expect the boundary 
layer to develop based upon the Reynolds number and the extent of suppression 
in the throat section. These predictions are consistent with the heat transfer 
results. 

In order to predict the boundary-layer development at higher chamber pres- 
sures, figure 21 was drawn. Within this figure, the distribution of K ax and 

Re@ are shown along the nozzle for various chamber pressures. K ax and Re© 
have several variables in common. Yet, K ax contains only freestream variables 
with the only boundary-layer variable, 0, contained in Reg. To understand the 
relaminarization phenomena, one must examine both these quantities. In Nash- 
Webber's study, Re© is linked to K by showing that as Re© increases, K cr 
increases. 

In figure 21, the momentum thickness Reynolds number was calculated using 
momentum thickness for a laminar boundary layer. Following each curve from 
left to right, the flow goes from the end of the combustor section to about 
0.5 in. past the throat. A dashed curve is drawn to indicate the position of 

the throat. As one goes from the upper curve to the lower curve, chamber pres- 

sure equals 656, 1004, 1700, 2100, and 2600 psia. This range corresponds to 
throat diameter Reynolds numbers of 5.6x1 0^ to 22x1 0 5 . 

Upon examining figure 21, one needs to bear in mind the approximate nature 
of these results. As mentioned previously, boundary-layer development depends 
on many factors such as the amount of acceleration, wall cooling, and rough- 
ness. The critical values of K ax and Re©, therefore, also depend upon these 
factors. Therefore, K ax cannot be used as an exact criteria, but, instead, as 
an indication of a trend or propensity towards a certain flow behavior. 

For the three upper curves, in the convergent section * *^ax is predicted 
to be quite high while Re© is relatively low. Therefore, laminar or transi- 
tional flow seems likely. For the two lower curves (2100 and 2600 psia), in 
the convergent section K ax borders on the critical value and Re© has risen 
above its critical value. Thus, transitional or turbulent flow could exist. 
Therefore, in the convergent section, the transition zone may extend up to 
2600 psia. For all of the curves, the acceleration parameter drops after the 
throat and the momentum thickness Reynolds number increases sharply. As the 
flow travels toward the exit, the acceleration parameter drops gradually and 
the momentum thickness Reynolds number remains at a fairly high value. There- 
fore, in the divergent nozzle, one would expect transitional or turbulent flow. 


CONCLUDING REMARKS 

Based upon the results of this analysis, the following recommendations 
are made. Because fully developed turbulent flow was not attained during this 
experiment, tests should be run on the nozzle at higher chamber pressures. 
Suggested run conditions are chamber pressures of 1500, 2000, 2500 and 
3000 psia. Fully developed turbulent flow will be indicated as the chamber 
pressure at which thrust or wall heat flux matches the turbulent prediction, or 
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the point at which thrust and wall heat flux remain at the same relative point 
between the laminar and turbulent prediction as chamber pressure increases. 

The latter condition is chosen because the prediction program may not be able 
to predict turbulent performance as accurately as necessary. Fully developed 
turbulent flow will also be indicated (as in figures 10 to 13) by a heat flux 
profile with a slope that is proportional to the chamber pressure to the 

0.8 power. Also, boundary-layer measurements should be taken to provide fur- 
ther data on the boundary-layer development. 

As an alternative to running the nozzle at chamber pressures high enough 
to produce turbulent flow in the convergent section, an analysis could be done 
to estimate the acceleration effects and the transition location for a specific 
chamber pressure. Then, the TDK program would be run while altering the tran- 
sition locations until the predicted heat flux distribution matches the experi- 
mental results. At that point, it would establish the amount of drag being 
produced which could be used in the nozzle performance database. 

Other options to produce fully developed turbulent flow in the nozzle at 
lower chamber pressures include reducing the convergence angle or increasing 
the combustor length to produce a stronger turbulent boundary layer in the con- 
vergent section. Also, roughening the wall of the convergent section would 
produce more turbulence. 


SUMMARY OF RESULTS 

The objective of this study was to determine the rocket engine chamber 
pressures at which boundary-layer flow in the 1030:1 nozzle is laminar, transi- 
tional and fully developed turbulent flow. Tests were performed on the nozzle 
at chamber pressures from 350-1000 psia. These conditions correspond to throat 
diameter Reynolds numbers from 3x10$ to 10xl0 3 * 5 . The propellants used were gas- 
eous oxygen and gaseous hydrogen. Thrust measurements and nozzle outer wall 
temperature measurements were taken during the 3-sec test runs. The hardware 
tested was a modified Rao contoured nozzle with a 25° half angle of conver- 
gence, a 1 in. throat diameter, a contraction area ratio of 4.223, and an exit 
area ratio of 1030. The copper throat section was water-cooled and expanded 
to an area ratio of 30:1. The carbon steel nozzle skirt was heat-sink. 

Results were based on a comparison of the experimental data and corre- 
sponding predictions from the Two-Dimensional Kinetics (TDK) nozzle analysis 
program. Predicted flow conditions were also used in conjunction with results 
from previous studies to estimate the state of the boundary layer at higher 
chamber pressures. The results are summarized as follows: 

1. Comparison of experimental and predicted wall heat transfer in the 
nozzle skirt indicated completely laminar flow at chamber pressures around 

360 psia, and transitional flow up to at least a chamber pressure of 1000 psia; 

2. Comparison of experimental and predicted vacuum thrust indicated lami- 
nar flow at a chamber pressure of 360 psia and laminar or transitional flow 

up to at least 747-psia chamber pressure. Thrust measurements at 1000-psia 
chamber pressure were not available; 

3. A least squares curve fit was applied to the wall heat flux data taken 

at an area ratio of 100 for various chamber pressures. Based upon an extra- 

polation of the curve, the transitional zone for the nozzle extends from a 
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throat diameter Reynolds number of 3 . 2x 1 0 5 to 22x1 0 5 - This corresponds to a 
chamber pressure range of 360 to 2600 psia; 

4. Nash-Webber 1 s lamlnarization map, which is used to predict boundary- 
layer development, indicated turbulent flow at a chamber pressure of 747 psia. 
Comparison with the experimental results suggests Nash-Webber 1 s criteria is not 
applicable to heat-sink or cooled hardware; 

5. Using predicted, flow conditions within the nozzle, distributions of 
the acceleration parameter K ax were calculated. These distributions 
indicated relaminarization effects in the convergent section which correlate 
with the heat transfer results found up to and including a chamber pressure of 
1000 psia. Relaminarization effects in the convergent nozzle were predicted 
up to a chamber pressure in the range of 2100 to 2600 psia, or throat diameter 
Reynolds number of 1 7 . 6x1 0^ to 22x10^. 

In conclusion, analysis of experimental data and predicted flow conditions 
indicates that for the NASA Lewis 1030:1 nozzle the laminar zone extends to a 
throat diameter Reynolds number of approximately 3.2xl0 5 . This corresponds to 
a chamber pressure of 360 psia. The transitional zone extends to a throat 
diameter Reynolds number of approximately 22xl0 5 , or a chamber pressure of 
2600 psia. Above this point, fully developed turbulent flow is predicted. 

These results are consistent with previous studies. 


REFERENCES 

1. Pavli, A.J., Kacynski , K.J., and Smith, T.A., "Experimental Thrust Per- 
formance of a High-Area-Ratio Rocket Nozzle," NASA TP-2720, 1987. 

2. Kacynski, K.J., Pavli, A.J., and Smith, T.A., "Experimental Evaluation of 
Heat Transfer on a 1030:1 Area Ratio Nozzle," NASA TP-2726, 1987. 

3. Schlichting, H., Boundary-Layer Theory . McGraw-Hill, New York, 1968. 

4. Kays, W.M., and Crawford, M.E., Convective Heat and Mass Transfer . 2nd 
edition, McGraw-Hill, New York, 1980. 

5. Moretti , P.M., and Kays, W.M., "Heat Transfer to a Turbulent Boundary 
Layer with Varying Free-Stream Velocity and Varying Surface Temperature 
-- An Experimental Study," International Journal of Heat and Mass Trans- 
fer . Vol. 8, Sept. 1965, pp. 1187-1202. 

6. Back, L.H., Cuffel , R.F., and Massier, P.F., "Laminarization of a Turbu- 
lent Boundary Layer in Nozzle Flow — Boundary Layer and Heat Transfer 
Measurements with Wall Cooling," ASME Paper 69-HT-56, Aug. 1969. 

7. Schoenman, L., and Block, P., "Laminar Boundary-Layer Heat Transfer in 
Low-Thrust Rocket Nozzles," Journal of Spacecraft and Rockets . Vol. 5, 

No. 9, Sept. 1968, pp. 1082-1089. 

8. Nash-Webber, J.L., "Wall Shear-Stress and Laminarization in Accelerated 
Turbulent Compressible Boundary Layers," Report 94, Gas Turbine Labora- 
tory, MIT, Apr. 1968. 


16 



9. Boldman, D.R., Schmidt, J.F., and Gallagher, A.K., "Laminarization of a 
Turbulent Boundary Layer as Observed from Heat Transfer and Boundary 
Layer Measurements in Conical Nozzles," NASA TN D-4788, 1968. 

10. "Liquid Rocket Engine Self-Cooled Combustion Chambers," NASA SP-8124, 

1977. 

11. Smith, T.A., Pavli, A.J., and Kacynski , K.J., "A Comparison of Theoretical 
and Experimental Thrust Performance of a 1030:1 Area Ratio Rocket Nozzle 
at a Chamber Pressure of 2413 kN/m 2 (350 psia)," NASA TP-2725, 1987. 

(Also AIAA Paper 87-2069). 

12. Nickerson, G.R., "The Rao Method Optimum Nozzle Contour Program," SEA 
Report 6/82/800.1, Software and Engineering Associates, Inc., Carson City, 
NV, June 1982. 

13. Evans, R.M., "Boundary Layer Integral Matrix Procedure - User's Manual," 
AEROTHERM-UM-75-64, Aerotherm Acurex Corp., Mountain View, CA, July 1975, 
NASA CR- 144046. 

14. Nickerson, G.R., Coats, D.E., and Dang, L.D., "Two-Dimensional Kinetic 
Reference Computer Program (TDK)," SN-63, Software and Engineering Associ- 
ates, Carson City, NV, Apr. 1985, NASA CR-178628. 

15. Keller, H.B., and Cebeci , T., "Accurate Numerical Methods for Boundary 

Layer Flows. II: Two-Dimensional Turbulent Flows," AIAA Journal, 

Vol . 10, No. 9, Sept. 1972, pp. 1193-1199: 

16. Cebeci, T., and Smith, A.M.O., Analysis of Turbulent Boundary Lavers , 
Academic Press, New York, 1974. 

17. Davidian, K.J., Dieck, R.H., and Chuang, I., "A Detailed Description of 
the Uncertainty Analysis of High Area Ratio Rocket Nozzle Tests at the 
NASA Lewis Research Center," NASA TM-1 00203, 1987. 

18. Praharaj, S.C., "Boundary Layer Simulator (BLIMPJ) Improvement," The 22nd 
JANNAF Combustion Meeting, Vol. 1, CPIA-PUBL-432, D.S. Eggleston, ed., 
Chemical Propulsion Information Agency, Johns Hopkins Univ., Laurel, MD, 
1985, pp. 477-493. (Avail. NTIS, AD-A165503). 


17 




L. 


CONVERGENT 


NOZZLE 


K. 


THROAT 


L 






TEST CAPSULE -i DIFFUSER 


EJECTORS 


SPRAY COOLER 


VERTICAL DRAIN LINE 


FIGURE 2. - ROCKET ENGINE TEST FACILITY (RETF) WITH CUTAWAY VIEWS OF THE TEST CAPSULE AND 
SPRAY COOLER. 
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FIGURE 3. - DIMENSIONS OF TEST HARDWARE (NOT TO SCALE) 



FIGURE 4. - ALTITUDE TEST CAPSULE - 1030:1 NOZZLE BEING INSTALLED IN THRUST STAND. 
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AXIAL DISTANCE FROM THROAT, IN. 

FIGURE 5. - EXPERIMENTAL AND PREDICTED WALL HEAT 
FLUX AT A CHAMBER PRESSURE OF 360 PSIA. 
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FIGURE 8. - EXPERIMENTAL AND PREDICTED WALL HEAT FLUX 
AT A CHAMBER PRESSURE OF 98A PSIA. . 
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HEAT FLUX 



CHAMBER PRESSURE, PSIA 

FIGURE 10. - EXPERIMENTAL WALL HEAT FLUX AS A FUNCTION OF 
CHAMBER PRESSURE AT AN AREA RATIO OF 100 (x = 7 IN.). 



CHAMBER PRESSURE, PSIA 


FIGURE 11. - EXPERIMENTAL WALL HEAT FLUX AS A FUNCTION OF 
CHAMBER PRESSURE AT AN AREA RATIO OF 300 (X = 1G.3 IN.). 
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HEAT FLUX 



FIGURE 12. - EXPERIMENTAL WALL HEAT FLUX AS A FUNCTION OF 
CHAMBER PRESSURE AT AN AREA RATIO OF 388 (X = 20.2 IN.). 



FIGURE 13. - EXPERIMENTAL WALL HEAT FLUX AS A FUNCTION OF 
CHAMBER PRESSURE AT AN AREA RATIO OF 800 (X = 39 IN.). 
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AXIAL DISTANCE FROM THROAT, IN. 

FIGURE 1G. - DISTRIBUTION OF K gx AT A CHAMBER PRES- 
SURE OF 360 PSIA. 
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AXIAL DISTANCE FROM THROAT, IN. 

FIGURE 17. - DISTRIBUTION OF K ax AT A CHAMBER PRES- 
SURE OF 656 PSIA. 
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AXIAL DISTANCE FROM THROAT, IN. 

FIGURE 18. - DISTRIBUTION OF K ax AT A CHAMBER PRES- 
SURE OF 748 PSIA. 
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FIGURE 19. - DISTRIBUTION OF K ax AT A CHAMBER PRES- 
SURE OF 984 PSIA. 
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AXIAL DISTANCE FROM THROAT, IN. 

FIGURE 20. - DISTRIBUTION OF K ax AT A CHAMBER PRES- 
SURE OF 1004 PSIA. 
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FIGURE 21. - PREDICTED VARIATION OF K ax AND MOMENTUM 
THICKNESS REYNOLDS NUMBER ALONG THRUST CHAMBER. 
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